Effective diffusivity of passive scalars in rotating flows 
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We use direct numerical simulations to compute turbulent transport coefficients for passive scalars 
in turbulent rotating flows. Effective diffusion coefficients in the directions parallel and perpendicular 
to the rotations axis are obtained by studying the diffusion of an imposed initial profile for the 
passive scalar, and calculated by measuring the scalar average concentration and average spatial 
flux as a function of time. The Rossby and Schmidt numbers are varied to quantify their effect on 
the effective diffusion. It is find that rotation reduces scalar diffusivity in the perpendicular direction. 
The perpendicular diffusion can be estimated from mixing length arguments using the characteristic 
velocities and lengths perpendicular to the rotation axis. Deviations are observed for small Schmidt 
numbers, for which turbulent transport decreases and molecular diffusion becomes more significant. 



I. INTRODUCTION 

The turbulent transport of passive scalars in 
anisotropic flows plays an important role in a wide variety 
of astrophysical and geophysical processes. The mixing 
in such flows is often studied and modeled by means of 
anisotropic transport coefficients. 

While the transport of passive scalars in isotropic tur- 
bulent flows has been received a lot of attention (see, e.g., 
[1-3]), less studies have considered the case of diffusion 
of the passive scalar in anisotropic flows. In the presence 
of anisotropy, the transport of the passive scalar is mod- 
ified and turbulent transport coefficients should become 
anisotropic. Of particular relevance for the atmosphere 
and the oceans (as well as for mixing in stellar convec- 
tive regions) are the cases or rotating and/or stratified 
turbulence. In the case of stratified flows, numerical sim- 
ulations [4| showed that the stratification acts to reduce 
the horizontal fluctuations of a scalar (i.e., it maintains 
the turbulent mixing in horizontal planes), while it sup- 
presses the vertical mixing and transport of the passive 
scalar. These results were later confirmed analytically in 
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Turbulent transport and diffusion of passive scalars 
in rotating turbulent flows has received less attention, 
although the anisotropy introduced by rotation is im- 
portant at the largest scales of the atmosphere and the 
oceans, as well as in structures in the atmosphere that 
develop strong local rotation, as in tornadoes [H, 0]- In 
the astrophysical context, the effects of rotation on tur- 
bulent diffusion are also important to understand the de- 
velopment latitudinal entropy gradients in rotating stars 
Furthermore, anisotropic passive scalar transport is 
believed to be associated with the partial depletion of 
lithium observed in the sun Q. 

The development of anisotropy in rotating flows dif- 
fers from that in stratified flows. While in the later case 
structures tend to be flat and layered, in rotating flows 
structures are elongated along the axis of rotation and 
the flow becomes quasi- two dimensional [To| . In rotating 



flows, turbulent transport coefficients for passive scalars 
were obtained from numerical simulations using the test 
field model in [8]. The results indicate that turbulent 
diffusion in the direction of rotation is reduced when the 
scales of the turbulence are comparable with those of the 
mean field. Recently, we considered scaling laws and in- 
termittency of the passive scalar in rotating flows [Tl| 
and found that the distribution of the passive scalar be- 
comes highly anisotropic, with scaling laws developing in 
the plane perpendicular to the rotation axis, and with 
scaling compatible with Kraichnan's model [12[ for the 
passive scalar in two dimensions. 

In a more general context, many authors have stud- 
ied turbulent diffusion in anisotropic flows in terms of 
particle dispersion. These studies are related with the 
transport of passive scalars, as the average concentration 
of a passive scalar is related to single-particle dispersion 
in a fluid. Theoretical and numerical studies [l3|, [lij (see 
also and simulations in [18j) indicate one parti- 

cle dispersion is suppressed along the direction of strati- 
fication in stably stratified turbulence with and without 
rotation. Moreover, results in [14] and [l7| show that the 
anisotropy introduced by stratification also reduces hor- 
izontal diffusion, although to a much lesser extent than 
vertical diffusion. These results are in good agreement 
with previous experimental results [19]. In the case of ro- 
tating flows without stratification, the same studies indi- 
cate that horizontal diffusion is highly reduced while ver- 
tical diffusion remains nearly unchanged. More recently, 
Lagrangian velocity autocorrelations (related with tur- 
bulent diffusion) were obtained in experiments of forced 
rotating turbulence [20|, and numerical simulations con- 
sidered mixing in rotating flows [21] . 

In this work, we study turbulent transport coefficients 
for the passive scalar in a turbulent rotating flow. Fol- 
lowing a procedure similar to the one used in [4[ for strat- 
ified flows, effective coefficients are obtained by studying 
the diffusion of an initial distribution of the scalar quan- 
tity in a turbulent flow, and calculated by measuring the 
average concentration and average flux of the scalar. Un- 
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like previous jobs (see, e.g., [8]), the coefficients obtained 
here are not scale dependent. Coefficients calculated in 
this fashion are useful when a computationally economi- 
cal description of the problem is required, as often is the 
case, e.g., in weather and climate modeling (see [6[). 

The velocity field and the passive scalar concentration 
are obtained from direct numerical simulations (DNS) of 
the Navier- Stokes equations in a rotating frame, together 
with the advection-diffusion equation for the scalar con- 
centration. All DNS are performed using 512 3 grid points 
in a regular periodic tridimensional grid. Three different 
rotation rates are considered, chosen to study turbulence 
with moderate Rossby number. Diffusion coefficients 
are calculated for the vertical (parallel to the rotation 
axis) and horizontal directions. A more detailed study 
of turbulent diffusion is then made varying the Schmidt 
number, in order to observe Schmidt and Peclet number 
effects. Finally, we present visualizations of the scalar 
which allow a more complete interpretation of the diffu- 
sive processes in the vertical and horizontal directions. 

We find that rotation dramatically reduces scalar diffu- 
sivity in the horizontal direction (in comparison with the 
isotropic and homogeneous case). This is in good agree- 
ment with theoretical arguments in [17] for single-particle 
dispersion, where it was found that vertical diffusion is 
twice larger than horizontal diffusion in the presence of 
pure rotation. Within error bars, our results are consis- 
tent with this finding. Moreover, we find that horizon- 
tal diffusion can be reasonably estimated using mixing 
length arguments as U±L± (where U± and L± are respec- 
tively the characteristic velocity and lengthscale in the 
direction perpendicular to the rotation axis). For small 
Schmidt and Peclet numbers, the turbulent diffusion de- 
creases as molecular diffusion becomes more important, 
while for large enough Peclet numbers the effective dif- 
fusion becomes independent of the Peclet number. 



II. NUMERICAL SIMULATIONS 

Data analyzed in the following section stems form DNS 
of the incompressible Navier- Stokes equations in a ro- 
tating frame for the velocity u, and of the advection- 
diffusion equation for the passive scalar given by 



d t u- 



• Vu = -2Q x u - Vp + iA7 2 u + f , (1) 



Vu = 0, 



d t + u • V<9 = kV 2 6. 



(2) 



(3) 



Here p is the pressure divided by the mass density (cho- 
sen uniform and constant in all simulations), v is the 
kinematic viscosity, and n is the scalar diffusivity. Also, 
f is an external force that drives the turbulence, and 
ft = Qz where Q is the rotation angular velocity. The 
mechanical forcing f is a superposition of Fourier modes 



with random phases, delta-correlated in time, injected in 
a narrow band of wavenumbers k G [1,2] (therefore we 
will consider the forcing wavenumber as kp ~ 1). 

Equations (pQ) , (|2]) and (|3)) are solved in a three dimen- 
sional domain of size 2ir with periodic boundary condi- 
tions using a parallel pseudospectral code [22|, [23| . The 
pressure is obtained by taking the divergence of Eq. (pQ), 
using the incompressibility condition (|2]), and solving the 
resulting Poisson equation. The equations are evolved in 
time using a second order Runge-Kutta method. The 
code uses the 2/3-rule for dealiasing, and as a result the 
maximum wavenumber is k max = N/3, where TV is the 
number of grid points in each direction (with N = 512 in 
all the runs). All simulations presented are well resolved, 
in the sense that the dissipation wavenumbers k v and k K 



are smaller than the maximum wavenumber k n 



at all 



times. 

The dimensionless numbers used to characterize the 
runs are the Reynolds, Peclet, and Rossby numbers, de- 
fined respectively as 
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where U is the r.m.s. velocity in the turbulent steady 
state, and L is the forcing scale of the flow defined as 
L = 2tt /kp- S c is the Schmidt number, defined as S c = 
v I k. In all simulations, U ~ 1, and the viscosity is v — 
6 x 10 -4 . The passive scalar diffusivity k is set equal to 
the viscosity for the main set of runs, although similar 
simulations but changing the value of S c (and therefore 
of P e ) were also performed. 

The detailed procedure followed in the numerical sim- 
ulations is as follows. We first conducted a simulation 
solving Eqs. (pQ) and (|2]) (Navier- Stokes without passive 
scalar and without rotation), starting from the fluid at 
rest (u = 0), and applying the external forcing f until 
reaching a turbulent steady state. This run was contin- 
ued for approximately 13 turnover times. The final state 
of this run was used as initial condition for the velocity 
field in multiple runs in which the external forcing f was 
kept the same but a passive scalar was injected. These 
runs can be grouped (namely in sets A, B, C, and D), 
with each set corresponding to several runs with the same 
rotation rate and Rossby number (see table [I]). 

Each run in each set corresponds to a simulation in 
which an initial Gaussian profile for the passive scalar 
was injected as follows: 



0(t = Q,Xi) = e-^-^ 2/a2 



(7) 



where i = 1 or 3 (i.e., the initial profile is a function of 
x\ = x or xs = z), fjL = 7r (the profile is centered in the 
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TABLE I: Parameters used in each set of runs. Q is the 
rotation rate, R Q is the Rossby number, R e is the Reynolds 
number, and U is the r.m.s. velocity in the turbulent steady 
state. 



Set 


n 


Ro 


Re Urms 


A 





oo 


1050 1 


B 


2 


0.04 


1050 1 


C 


4 


0.02 


1050 1 


D 


8 


0.01 


1050 1 


TABLE II: Parameters 


used in one 


of the subsets of runs. 


In the name of each run, the subindex x or z indicates the 


dependence with the coordinate of the initial passive scalar 


Gaussian 


profile, and the subindex with the number is 1/S C > 


For each 


run, v is the kinematic viscosity, n is the molecular 


diffusivity, P e is the Peclet number, and S c is the Schmidt 


number. 








Run 


V 


ft 


Pe Sc 


A x i/2 


6 x 10" 4 


3 x 10~ 4 


2100 2 


A xl 


6 x 10" 4 


6 x 10" 4 


1050 1 


A x2 


6 x 10" 4 


1.2 x 10" 3 


520 0.5 




6 x 10" 4 


2.4 x 10" 3 


260 0.25 


A x8 


6 x 10" 4 


4.8 x 10" 3 


130 0.12 


A X IQ 


6 x 10" 4 


9.6 x 10" 3 


66 0.06 


A X 32 


6 x 10" 4 


1.92 x 10" 2 


33 0.03 



scalar, we consider quantities averaged over the two di- 
rections perpendicular to the initial dependence of the 
Gaussian profile. In particular, we consider the average 
passive scalar concentration and the spatial flux Oiti, 
where i = 1 or 3 depending on the initial dependence of 
the Gaussian profile, and the averages are done over the 
two remaining Cartesian coordinates. Note the spatial 
flux Qui represents the amount of passive scalar trans- 
ported in the i-direction per unit of time by the fluc- 
tuating (or turbulent) velocity. Since there is no mean 
flow in the simulations (we use delta-correlated in time 
random- forcing), ui is the fluctuating velocity. Then, the 
effective turbulent diffusion coefficient is Q 



Di(xi) 



(8) 



This coefficient corresponds to how much passive scalar is 
transported by the fluctuating velocity, per unit of varia- 
tion of 6 with respect to X{. As already mentioned, i = 1 
will stand for horizontal diffusion, while i=3 will stand 
for vertical diffusion. 



III. NUMERICAL RESULTS 



A. Isotropic turbulence 



middle of the box, with the box of length 27r), and a = 1. 
This allows us to study the diffusion of the initial profile 
in the direction parallel to rotation (z, or vertical) and 
in the direction perpendicular (x, or horizontal). In a 
few runs it was verified that the diffusion in the x and y 
directions was the same (as rotating turbulence tends to 
be axisymmetry ) . The runs in each set are labeled with a 
subindex indicating the dependence of the initial profile 
(e.g., runs in group A are labeled A x or A z depending on 
the initial Gaussian profile used). 

Finally, to study the effects of varying the Schmidt 
number, several simulations were done increasing the 
molecular diffusivity from ft = v (all the simulations de- 
scribed above) to ft = 32v (increasing 1/S C by factors of 
two in each run). Simulations with ft = v/2 were also per- 
formed. To measure the effective diffusivity all simula- 
tions in each group (with subscript x or z) were repeated 
with different values of ft. This results in twelve more 
simulations in each group, with S c = 2 to S c = 1/32. 
To differentiate each one in each set, a subindex equal to 
1/S C is added after the subindex x or z. The list of all 
runs in set A is shown in table Ull indicating the value 
of ft, and the corresponding P e and S c numbers. The 
runs in the other sets are labeled following the same rules, 
and have the same parameters as the ones shown in Ta- 
ble HI As an example, runs A X 4 and C X 4 have the same 
values of v, ft, P ei and S c = 1/4, with an initial Gaussian 
profile for the passive scalar in x. The two runs differ 
though in the Rossby number and rotation rate. 

To characterize the turbulent diffusion of the passive 



In the absence of rotation (runs in set A), we expect 
the diffusion coefficients to be isotropic (i.e., the horizon- 
tal and vertical diffusions to be the same within error 
bars). Figure Q] shows the horizontal average profile of 
the passive scalar 0(x), the horizontal flux 6u x (x), and 
the horizontal diffusion D x (x) at five different times in 
run A x \. As time evolves, the mean profile 6 flattens 
and widens. The flux is roughly antisymmetric, and is 
positive for x > tt and negative for x < tt. This can be 
expected as there is an excess of 6 at x = tt at t = which 
turbulent mixing should diffuse, by transporting this ex- 
cess towards x = and x = 2tt. Horizontal diffusion 
increases to a maximum value around t ~ 1.5; after this 
time it fluctuates around an average value. Large fluc- 
tuations in D x (x) near the center of the box are due to 
the fact that, by definition, turbulent diffusion diverges 
at that point. 

Passive scalar concentration, vertical flux and vertical 
turbulent diffusion were also calculated for run A z ± (i.e., 
the same run but with an initial Gaussian profile in z). As 
expected, the same results as in run A x \ were obtained. 

B. Effect of rotation 

1. Horizontal transport 

Figure [2] shows the mean horizontal concentration 0(x), 
the horizontal flux 0u x (x), and the horizontal diffusion 
D x {x) for different times in run C x \ (fi = 4). In this case 
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FIG. 1: (a) Average horizontal concentration # in run A x i, 
at times £ = 0, 0.5, 1 and 1.5, denoted by solid, dotted, dash, 
and dash-dotted lines respectively, (b) Horizontal flux at the 
same times, (c) Horizontal turbulent diffusion at the same 
times. 



FIG. 2: (a) Averaged vertical concentration in run C x i, at 
times t = 0,0.25,0.5,0.75, and 1. denoted by solid, dotted, 
dash, dash-dotted and dash-triple-dotted lines respectively, 
(b) Horizontal flux at the same times, (c) Horizontal turbu- 
lent diffusion at the same times. 



diffusion grows faster at earlier times, but then saturates 
at a lower value. Also note that the average profile and 
the average flux become asymmetric: there is an excess 
of 6 for x < 7r, and the flux (in absolute value) is larger 
for x < 7r than for x > tt. This asymmetry is associated 
with the Coriolis force and has been already observed in 
Q. As will be shown later, it results from a rotation of 
the initial passive scalar profile. 

From D x (x) we can obtain a global turbulent diffusion 
coefficient D x by computing its average over all values 
of x. Figure [3] shows the resulting horizontal turbulent 
diffusion as a function of time for runs A x \ and C x \. 



In both runs D x grows from an initially small value to 
a saturation value between t « 1 to « 1.5. Later D x 
displays large fluctuations, and if the runs are continued 
further it decreases as the profile of the passive scalar is 
finally diffused and homogeneity is recovered. 

Although at early times D x grows faster in the run 
with rotation (run C x i), at later times D x saturates at 
a lower value indicating horizontal diffusion decreases in 
the presence of rotation (in comparison with the isotropic 
and homogeneous case). This is in good agreement 
with predicted reductions of the horizontal passive scalar 
transport in the presence of rotation [13, EH, and also 
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FIG. 3: (a) Horizontal turbulent diffusion as a function of 
time for runs A x i (solid) and C x i (dotted), (b) Detail in log- 
log coordinates of the horizontal turbulent diffusion for run 
A x i. (c) Same for run Cxi • 



with previous measurements of turbulent transport coef- 
ficients in numerical simulations using a different method 
(the test field model, see [8[). 

It is interesting that at intermediate times the growth 
of D x can be roughly explained by a simple phenomeno- 
logical model. If we think that turbulent eddies dominate 
the transport over molecular diffusion, and we think of 
the turbulent flow as a superposition of eddies at differ- 
ent scales £, then we can imagine that the smallest ed- 
dies (much smaller than the integral scale eddies), with 
faster turnover times rg, are the first to start mixing the 
passive scalar. As time evolves, larger and larger eddies 



FIG. 4: Horizontal turbulent diffusion as a function of the 
inverse Rossby number, obtained from simulations A x i, B x i, 
Cxi , and D x i (all runs with S c = 1)- Error bars correspond 
to the mean standard deviation. 

come into play, as the eddies at larger scales are able to 
complete a turnover. For eddies in the inert ial range, the 
eddy turnover time can be estimated as 

n ~ £/ui (9) 

In isotropic and homogeneous turbulence, the inertial 
range scaling for the velocity field (in the absence of inter- 
mittency corrections) can be written as ug ~ i 1 / 3 . Then, 
the eddy turnover time is 

n~£ 2/3 . (io) 

From mixing length arguments, the turbulent diffusion 
at late times is D x ~ LU. At early times, if only the 
eddies with turnover time smaller than the actual time 
contribute to the mixing, then from Eqs. (|9)) and (fTU|) we 

get 

D x ~rl (11) 

In the rotating case the inertial range scaling for the 
velocity field is modified. Rotation sets a preferen- 
tial direction for energy transfer, resulting in a quasi- 
bidimensionalization of the flow. Assuming E(k±) ~ kj 2 
in the inertial range (see, e.g., [24], [25|), it follows that 

1 /2 

U£ ± ~ l_[ . If at late times we can assume the horizon- 
tal diffusion D x ~ L±U± (with L±_ and U± respectively 
the characteristic length and velocity in the perpendic- 
ular direction), and using the same arguments as in the 
isotropic and homogeneous case, then 

D x ~ T f m (12) 

Assuming t ~ t$ (i.e., that at a given time only the eddies 
that were able to do a full turn contribute to the mix- 
ing) , these scalings are approximately fulfilled by the runs 
right before D x reaches saturation, and are indicated as 
a reference in Fig. [3j 

The reduction of the horizontal turbulent diffusion ob- 
served in Figs. [TJ [2j and [3] can be further confirmed by 
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FIG. 6: Vertical turbulent diffusion as a function of the inverse 
Rossby number, obtained from simulations A z ±, B z ±, C z i, and 
D z i (all runs with S c = 1)- Error bars correspond to the mean 
standard deviation. 
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FIG. 7: Turbulent diffusion D x as a function of the Peclet 
number for runs J\ x \ to A X 32 (H = 0). Note that as the flow 
in these runs is approximately isotropic, D x w D z . 



FIG. 5: (a) Averaged passive scalar profile 6(z) in run C z i, 
at times t = 0, 0.25, 0.5, 0.75, and 1, denoted by solid, dotted, 
dash, dash-dotted and dashed-triple-dotted lines respectively, 
(b) Mean vertical flux as a function of z at different times in 
the same run. (c) Same for the vertical turbulent diffusion. 



computing the time average of D x during the steady state 
of all runs with S c = 1 . Figure |4] shows the mean value 
of the horizontal diffusion as a function of the Rossby 
number for runs A x i, B x i, C^i, and D x \ (i.e., simula- 
tions with ft = 0,2,4 and 8, and with initial Gaussian 
profile of the scalar in x). The error bars correspond to 
the mean standard deviation. It is clear from the data 
that D x decreases monotonically as 1/R Q increases. 



2. Vertical diffusion 



Figure [5] shows the mean vertical passive scalar concen- 
tration 0(z), the mean vertical flux 0v z (z), and the ver- 
tical diffusion D z (z) at different times in run C z ±. Note 
the profiles here are more similar to those obtained in 
the isotropic and homogeneous case: 0{z) and 0v z {z) are 
respectively symmetric and antisymmetric with respect 

tO Z = 7T. 

As in the case of horizontal diffusion, we can compute 
the mean vertical diffusion coefficient in the steady state 
of runs A z \, B z \^ C z \, and D z \. This results in a depen- 
dence of D z with the Rossby number, and is shown in 
Fig. [6] While horizontal diffusion is strongly dependent 
on the Rossby number, the vertical diffusion seems to be 
(within error bars) independent of its value. 
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FIG. 8: Horizontal turbulent diffusion as a function of the 
Peclet number for runs C x \ to C x ?>2 (H = 4). 
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FIG. 10: (a) Horizontal energy spectrum (solid) and passive 
scalar variance spectrum (dash-dotted) at late times in run 
Cxi- (b) Same for run C X 32- 



FIG. 9: Vertical turbulent diffusion as a function of the Peclet 
number for runs C z \ to C z ?>2 (Q = 4). 



C. Effect of S c and P e numbers 

A large number of runs were performed to study the 
effect of the Schmidt number in turbulent diffusion. As 
the molecular diffusivity n is increased, we can expect 
turbulent diffusion to decrease until molecular diffusion 
dominates. Figure shows the turbulent diffusion D x as 
a function of P e in runs with Q = (i.e., estimating D x 
from all the runs in set A). As the flow is approximately 
isotropic, D x « D z . 

Effects associated with the Peclet and Schmidt num- 
bers can be observed for small P ei as D x starts to de- 
crease for P e < 100. However, for a wide range of values 
of the Peclet number D x remains approximately con- 
stant. This is in good agreement with theoretical ex- 
pectations that for small enough v and ft, the effective 
turbulent Schmidt number should be of order one (see 
(26^ for similar arguments in the case of the turbulent 
magnetic Prandtl number). 

A similar behaviour is observed in the rotating case, as 
can be seen from Figs. [8] and [9] In these figures, the tur- 
bulent transport coefficients are computed from the runs 
in set C (Q. = 4, R Q = 0.02). As in the isotropic case, for 



sufficiently large P e turbulent diffusion remains approx- 
imately constant and the effective Schmidt number is of 
order one. P e number effects are observed earlier for ver- 
tical diffusion at small P e (indeed, for P e < 500 a change 
in the behaviour of D z can already be seen), although 
it seems that much smaller P e numbers are required to 
see a significant decrease in the horizontal diffusion. This 
could be associated with the fact that in the presence of 
rotation the turbulent diffusivity is already reduced, and 
that therefore larger molecular diffusivities are needed to 
decrease it further. Also, note that in the runs with large 
enough P e and small enough R Q1 the vertical turbulent 
diffusion reaches a value that is approximately twice the 
horizontal, in good agreement with predictions in [l4l[l7j. 

The change in the turbulent diffusivity as the P e num- 
ber is decreased can be also understood from the spectra 
of energy and passive scalar variance (see Fig. [T0|) . For 
Sc = 1, at late times the energy and the passive scalar 
variance display an inertial range in a similar range of 
wave numbers. As a result, it can be expected that the 
turbulent eddies will mix the passive scalar concentration 
at all scales. For small S c and P e , the passive scalar can- 
not develop an inertial range. As a result, eddies smaller 
than the dissipation scale for the passive scalar (i.e., ed- 
dies in the inertial range of the velocity field) are unable 
to mix the passive scalar and thus the turbulent diffusion 
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FIG. 11: Passive scalar concentration in a horizontal plane in 
run Cxi, at times t = 0, 1, 1.5, and 2.5 from left to right and 
from top to bottom. 
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FIG. 12: Passive scalar concentration in a horizontal plane at 
t = 1 in runs A x i (left, no rotation) and D x i (right, Q = 8). 



decreases. 



D. Turbulent structures and diffusion 

From the results presented above it is clear that rota- 
tion plays an important role in the diffusion of the passive 
scalar, modifying its mixing with respect to the isotropic 
and homogeneous case. Figure [11] shows a cut in a hor- 
izontal plane of the passive scalar concentration in run 
Cxi at different times. The initial Gaussian profile not 
only diffuses in time, but also bends and rotates. The 
bending of the passive scalar concentration was observed 
before in [8| and explained as an effect of the Coriolis 
force. In our runs, the passive scalar at t = is con- 
centrated in a narrow band around x = 7r. The average 
flux is thus towards positive values of x for x > 7r, and 
towards negative values of x for x < ir (i.e., in the direc- 
tion of -V<9, see Fig. [2]). The Coriolis force in Eq. (pQ) 
is —2Vtz x u and therefore on the average this force cre- 



FIG. 13: Passive scalar concentration in a vertical plane in 
run Czi, at times t = 0, 1, 1.5, and 2.5 from left to right and 
from top to bottom. 



ates a drift of the flux towards positive values of y in the 
x > 7r region, and towards negative values of y for x < tt 
jS|. This explains the bending of the initial profile we 
observe of the runs with rotation, that is not observed in 
the runs without rotation (see Fig. [12] for a comparison). 

Diffusion in the parallel direction is of a different na- 
ture, and more strongly dependent of the structures 
that emerge in rotating turbulent flows. Rapidly ro- 
tating flows are characterized by columnar structures in 
the velocity field and vorticity, associated with a quasi- 
bidimensionalization of the flow. The mechanism un- 
derlying the transfer of energy towards two dimensional 
modes and responsible for the formation of these columns 
seems to be associated with wave resonances in the 
energy-exchanging triadic interactions (27J . Two-point 
closures of turbulence, such as the E ddy Damped Quasi- 
Normal Markovian closure (see, e.g., [24[) successfully ex- 
plain the emergence of columns with the same principle. 
However, there are alternative theories that consider the 
formation of columns as the result of a relative concentra- 
tion of kinetic energy in cylindrical structures resulting 
from the conservation of linear and angular momentum 

Columnar structures have been reported in many nu- 
merical simulations of turbulent flows (see, e.g., [29]). As 
these columns live for long times and move across the do- 
main, they play an important role in the mixing of the 
passive scalar. Figure [13] shows a cut in a vertical plane 
of the passive scalar concentration at different times in 
run C z \. Note that diffusion is different from the one ob- 
served in horizontal planes in the same run (Fig.[TT]h and 
from the one observed in the isotropic and homogeneous 
case (Fig. [12] (a)). The passive scalar is diffused from its 
initial profile in vertical stripes, that are stretched fur- 
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ther (thus increasing the mixing) as time evolves. This 
stripes are created by updrafts or downdrafts inside the 
columns. As these columns go through the region with 
large concentration of the passive scalar, the updrafts or 
downdrafts mix the passive scalar with the regions im- 
mediately above or below. 

IV. CONCLUSIONS 

We used 56 direct numerical simulations with regular 
spatial resolution of 512 3 grid points to measure turbu- 
lent diffusion in directions parallel and perpendicular to 
the rotation axis, in turbulent flows at different Rossby 
and Schmidt numbers. The effective coefficients were ob- 
tained by studying the diffusion of an initial concentra- 
tion of the passive scalar and calculated by measuring its 
average concentration and average spatial flux. 

The effect of rotation in turbulent diffusion is the 
op pos ite of that found in the presence of stratification 
[4j, LL8| : while in the latter case stratification reduces ver- 
tical diffusion with respect to horizontal diffusion, rota- 
tion dramatically reduces scalar diffusivity in the hori- 
zontal direction. In our simulations, vertical diffusion re- 
mains of the same order as in isotropic and homogeneous 



turbulence, although for a different reason: in rotating 
flows, it was found that diffusion in the vertical direc- 
tion is strongly dominated by updrafts and downdrafts 
in columnar structures in the velocity field. 

Within error bars, and for small enough Rossby num- 
ber and large enough Schmidt number, our results 
are consistent with theoretical results based on single- 
particle dispersion that predict that vertical diffusion is 
twice larger than horizontal diffusion in the presence of 
pure rotation [3, E3- For small Schmidt and Peclet 
numbers the turbulent diffusion decreases as molecu- 
lar diffusion becomes more important, while for large 
enough Peclet numbers the effective diffusion becomes 
independent of the Peclet number, resulting in a turbu- 
lent Schmidt number of order one. 
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